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A review is presented of the current status of experimental searches for physics beyond the standard model in 
rare r decays. 



1. Introduction 

The availability of large samples of r leptons al- 
lows for searches of very rare (Standard Model forbid- 
den) r decay modes, which provide a sensitive probe 
for physics beyond the Standard Model (BSM). In re- 
cent years, very large samples of r leptons have been 
made available by the high luminosity running of the 
two asymmetric 5-factories PEP — II and KEK £?, 
which, in light of the comparable size of the e + e~ — > bb 
and e + e~ — > r + r~ cross sections (cr r + r - ~ 0.9 nb, at 
y/s ~ Mr(45)), can rightfully be regarded as r fac- 
tories as well. In fact, the luminosity so far recorded 
by the experiments BABAR and Belle running respec- 
tively at PEP -II and KEK B (C BaBar ~ 470 fb" 1 , 
and C Belie ~ 710 fb -1 ) entails that the total sample 
of r pairs produced in the two experiments is already 
currently in excess of 10 9 . 

Differently from the quark sector, the lepton world 
is well known to be characterized by an apparent re- 
markable conservation of flavor. An exact conserva- 
tion of Lepton Flavor (LF) is naturally implemented 
in the Standard Model through the assumption that 
neutrinos are mass degenerate (actually: massless) 
fermions. Since the observation of neutrino oscilla- 
tions, we know that the latter assumption is actu- 
ally wrong, and that, with it, we also have to give 
up with the strict conservation of LF. However, the 
minimal extension of the SM that amounts to simply 
allowing for massive and nondegenerate neutrinos 1 , 




Figure 1: Dominant SM diagram contributing to the ra- 
diative LFV decay r — > fij. 



1 I will henceforth refer to such "naive" minimal extension of 




Figure 2: One loop SM diagram contributing to the non- 
radiative LFV decay r~ — > ji~l + l~ (from [1]). 



due to the remarkable effectiveness of the ensuing lep- 
tonic GIM mechanism which in turn stems from the 
smallness of the neutrino masses, predicts extremely 
low rates for all LF violating (LFV) phenomena in 
the charged lepton sector, notwithstanding the cur- 
rent evidence in favor of large leptonic mixing. Using 
the largest value in the currently favored experimen- 
tal range [J for Am| 2 = |ra§ - m\\ ~ 3 10" 3 eV 2 , 
where m^, i = 1,2,3 are the three neutrino mass 
eigenvalues, and allowing for maximal mixing, the 
SM expected rate [2j |3] for the radiative LFV decay 
r — > fij (see Fig. [I] for the dominant SM diagram) 
amounts to a surely unobservable B(r~ — > /jL~j) = 

ilt (^f ) 2 sin 2 20 23 B(r~ - ^V,v T ) ~ 10" 54 . 

It is worthwhile to remark that the SM suppres- 
sion of LFV processes in the charged sector, though 
a general feature of the model, doesn't achieve to 
the same degree, and in particular to the extreme 
one mentioned above, for all processes. As noted by 
Pham [4], LFV r decays of the kind of r~ — > fi~l + l~ 
or r~ — > fi~ p° (I = e, /i), due to the presence of in- 
frared divergences in the contributing diagram shown 
in Fig. |2j get in fact a milder, logarithmic, suppres- 
sion factor in place of the striking factor (Am^/M^) 
suppressing the class of radiative LFV processes r — > 
Ij. The SM prediction for the rate of such nonradia- 
tive LFV processes turns in fact to be proportional 



the Standard Model as "SM" , without further comments upon 
it actually being an extension of the Standard Model proper. 
It amounts to simply extending the Standard Model by the 
addition of an appropriate neutrino mass matrix in a gauge in- 
variant manner, while keeping the minimal Higgs scheme, with 
only one Higgs doublet. 
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to J2k=2 U T kU* k In Taking into account the cur- 
rent experimental knowledge on neutrino mixing, this 
translates to expected branching ratios many orders of 
magnitude larger than those expected for for r — > £7 
(Pham [4 estimates 6 > 10 -14 ), but in any case far 
below the current experimental reach. 

The above considerations allow us to draw two im- 
portant conclusions. The first one is that, although 
the discovery that neutrinos are massive entails that 
LF numbers are no longer conserved in the SM, we 
can safely state that, as long as the SM represents the 
correct physical description, no LFV process in the 
charged sector (CLFV) should ever be observed. This 
means that observation of CLFV would undoubtedly 
be a signature of BSM physics. The second conclusion 
is that it is worthwhile and important to search for 
CLFV in all possible experimentally searchable pro- 
cesses (e.g. radiative and nonradiative ones), since: 
(a) as shown to be the case in the SM, also in most 
BSM physics scenarios different CLFV processes could 
be characterized by widely different rates, and (b) the 
eventual observation and measurement of LFV in dif- 
ferent processes can be generally expected to be sensi- 
tive to different sectors of the new physics parameter 
space and better constraining the choice of a new sce- 
nario. 

In fact, many scenarios of BSM physics predict rates 
for LFV r decays that are not too far, or even within, 
the current experimental reach. Without pretending 
to be exhaustive, table [I] shows some examples of sce- 
narios that have been theoretically investigated, with 
the corresponding order of magnitude expectations for 
the branching ratio of radiative and nonradiative LFV 
r decays. 

Table I Examples of BSM scenarios predicting rates of 
LFV r decays not too far, or even within, the current ex- 
perimental reach. Order of magnitude expectations for the 
two broad classes of radiative (e.g. r — > Z7) and nonradia- 
tive (e.g. r — > 3/) LFV r decays are displayed in the last 
two columns. 
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2. Searching for LFV r decays at e + e~ 
colliders 

The BABAR [13 and Belle [14 detectors are re- 
markably similar, with the major difference being 



SicgTEjj decay 




/ SM x decay 
* (1- or 3- prong) 
(lor 2 v) 

Figure 3: Schematic drawing of the production of r pairs 
in e + e~ collisions, displayed in the CM frame. The figure 
also illustrates the general strategy employed in searches 
for LFV r decays at e + e~ colliders: one r is reconstructed 
in a SM allowed ("Tag") decay, characterized by a 1- or 3- 
prong + 1 or 2 v final state, while the search for a possible 
LFV (neutrinoless!) decay is performed on the recoiling 
("Signal") r. 



in the technology used to identify charged particles: 
Belle uses a threshold Cerenkov detector together 
with time-of- flight and tracker dE / dx, whereas BABAR 
mainly relies on a ring-imaging Cerenkov detector 
augmented by dE/dx in the trackers. 

In addition to the statistics issue already discussed 
in the introduction, e + e~ colliders running at, or 
nearby, the T(45) resonance provide a very clean en- 
vironment for the search of rare (and, in particular, 
LFV) r decays, r leptons are exclusively 2 produced 
in pairs through the QED process e + e~ — > r + r~, and 
thereby fly back-to-back in the e + e~ Center of Mass 
(CM) frame (Fig. |3|. At the asymmetric 5-factories 
PEP — II and KEK B the events are boosted in 
the laboratory (Lab) frame with a /?7 ~ 0.56 for the 
BABAR experiment and ~ 0.43 for Belle. The boost of 
each of the two r's in the CM frame, combined with 
such a CM-to-Lab boost, confer to the final state of 
each r decay a "jet-like" shape which allows not only 
to powerfully reject the BB background, but also to 
easily separate the final state products of one r from 
those of the other. 

The general strategy employed in searches for LFV 
r decays at e + e~ colliders consists in searching for 
e + e~ — > r + r~ events where one r is reconstructed in 
a SM allowed decay (henceforth referred to as "Tag" 



2 Apart from possible initial or final state radiation 7('s). 
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Figure 4: Scatter plot of the reconstructed variables 
Mlfvd and AE for MC simulated e + e~ — > t + t~ events 
where one of the r leptons decays to the LFV final state 
fij (from the BABAR experiment). For the exact meaning 
of the variable plotted along the vertical axis (Mec), see 
the text and footnote [3] The dark green ellipse drawn on 
top of the distribution in correspondence of the constraints 
expressed by Eqs. [I] and [2] illustrates in a qualitative fash- 
ion the relative size of a typical 2 — 3 standard deviations 
wide signal box employed in the final event selection. 
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Figure 5: Scatter plot of the reconstructed variables p m iss 
and ra^iss f° r data (dots) and signal MC (shaded boxes) 
events in a search for the LFV decay r — > /jj by the Belle 
collaboration [18] . The selection criterion employed con- 
sists in requiring that the event lies between the two lines 
in the figure. 



decay), characterized by a 1- or 3-prong + 1 or 2 v 
final state, while the search for a possible LFV decay 
(LFVD) is performed on the recoiling ("Signal") r. A 
noteworthy feature of all LFV r decays searched for 
by BABAR and Belle is that such decays do not bear 
any undetectable particle (specifically: any neutrino) 
in their final state, and can therefore be completely re- 
constructed. In other words, the observed missing mo- 
mentum of the reconstructed complete event is wholly 
attributable to the "Tag" decay. This circumstance 
provides very powerful tools for background rejection, 
which amount to requesting that the reconstructed in- 
variant mass and CM energy for the candidate LFVD 
( "Signal" ) side of the event coincide with the nominal 
mass of the r lepton and with half the CM energy of 
the e + e~ collision, respectively: 

M LFV D = m T (1) 

AE = E£ F M VD -^=0 (2) 

An example of the distribution of simulated signal 
(r — > fij) MC events in a scatter plot of these two 
kinematical variables is shown in Fig. [4j which dis- 
plays the accumulation of signal events in correspon- 
dence of the two constraints given by Eqs. [I] and [2j as 
well as the smearing brought about by resolution and 
radiative effects, respectively. The same figure dis- 
plays also in a qualitative fashion a two dimensional 



signal box which is used to separate the signal from 
the SM backgrounds. 

Exploiting at best the available constraints can sig- 
nificantly boost the sensitivity of the search. An ex- 
ample of this important remark is offered by the strat- 
egy employed by the BABAR experiment since the pub- 
lication of their search for r — > fij decays [15 . The 
choice of evaluating Mlfvd by imposing a beam en- 
ergy constraint 3 , obtained by fixing the "Signal" r de- 
cay vertex - or equivalently the emission point of the 
photon in the final state - in correspondence of the 
point of closest approach of the signal \i trajectory 
to the beam axis, brought a dramatical improvement 
in terms of resolution on such a kinematical variable 4 
with respect to a more straightforward reconstruction 
of the same quantity as an unconstrained invariant 
mass of the LFVD decay products. Typical values for 
the resolution on the other kinematical quantity AE 
range around ~ 50 MeV. 

Other ingredients crucial to an effective reduction 
of backgrounds consist in a clever use of the powerful 



3 Whence the use of the symbol Mec in the scatter plot of 

Fig J 

4 The resolution improvement on Mlfvd due to the beam 
energy constraint technique (Mec) is greater of a factor two: 
a (M E c) ~ 9 MeV, while cr(M inv ) ~ 20 MeV, where M inv de- 
notes the straightforward evaluation of the kinematical quan- 
tity Mlfvd as an unconstrained invariant mass of the LFVD 
decay products. 
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Figure 6: Schematic representation of a doubly radiative, 
e + e~ — > /i + /x~77, fake missing momentum event. 



particle identification capabilities of the two experi- 
ments and, whether possible, in the exploitation of 
additional constraints on the event. An example of 
the latter strategy is represented by the requirement 
that the reconstructed event missing mass (m m i ss ), 
when the "Tag" decay is chosen so as to have a single 
neutrino in the final state, be compatible with zero. 
Such a constraint proved to be very effective in Belle's 
searches for the LFV decay r — > /ry [18 , as illustrated 
in Fig. [5J where the selection criterion on m^ iss em- 
ployed by the Belle collaboration is seen to be imple- 
mented in a Pmisg-dependent fashion. 

In addition, a rewarding strategy has in some cases 
turned out to consist in trying to identify, and then 
cleverly fight, specific, particularly nasty, backgrounds 
which could otherwise survive all other selection cri- 
teria. An example of this situation is represented by a 
fraction of the residual /i/i7(s) background in searches 
for the r — > /ij decay, namely that consisting of dou- 
bly (or more) radiative events where the flight path 
of one of the photons lies sufficiently close to the tra- 
jectory of the charged track in the "Tag" side of the 
event. In these occurrences, schematically illustrated 
in Fig. |6j such a photon may fail to be reconstructed 
as a neutral cluster in the electromagnetic calorime- 
ter and its energy deposit in the calorimeter could 
eventually get associated to the nearby charged track. 
When the photon is sufficiently energetic, such occur- 
rences would mimic a substantial missing momentum. 
Concurrently, the "Tag" side charged track would er- 
roneously fail to be identified as a muon according to 
basically any of the standard /i identification criteria 
because of the large energy deposit in the calorimeter 
erroneously associated to it. 



As a general rule, both BABAR and Belle optimize 5 
their selection strategy, and the particular values of 
the associated cuts, by using samples of fully simu- 
lated and reconstructed signal and background MC 
events. MC samples, possibly checked with control 
samples or with data events sufficiently far from the 
signal region (sidebands), are also used to model the 
background shapes of relevant kinematical variables 
for the purpose of estimating and subtracting the fi- 
nal background contribution, whose normalization is 
obtained from sidebands data. The MC simulation 
of the signal is also used to determine the signal effi- 
ciency (e), which typically lies between 2% and 10%, 
depending on the channel under study. To minimize 
possible biases, both experiments adopt a blind anal- 
ysis approach, which consists in excluding from con- 
sideration, up to the moment when the full selection 
and all systematic studies are completely finalized, all 
events in the data within a suitably shaped and sized 
region in the Mlfvd and AE plane around the signal 
signature values given by Eqs.[l]and[2| 

In order to possibly boost the sensitivity of the 
search, the two experiments have adopted a spectrum 
of different approaches to finally discriminate a pos- 
sible signal from the residual background and obtain 
an estimate of the former's size: these range from sim- 
ple cut-and-count to likelihood techniques, and in one 
case a Neural- Network selection has been adopted. If 
the estimated background (AT bkd ) is compatible with 
the observed number of events, a 90% confidence level 
(CL) upper limit on the number of signal events is es- 
timated 6 (Nq l ). Based on this quantity and on the 
estimate of the signal efficiency e, a 90% CL upper 
limit on the branching ratio is obtained as: 



UL 
90 



N UL 
iV 90 



2Co~ r + r - e 



(3) 



3. Experimental results 

Both BABAR and Belle have searched for LFV in 
many different classes of r decays. In the following 
I'll summarize the results obtained by the two exper- 
iments, grouping them according to the decay final 
state searched for and giving particular emphasis only 
to the most recent ones. 



5 The analyses are optimized to give the best "expected up- 
per limit". 

6 Also in this respect, the two experiments have resorted to a 
number of different techniques 16 17 . It should be noted that 
the experimental upper limits quoted by both the experiments 
are in any case of a frequentistic nature. Negative fluctuations 
are generally admitted (and occur!) and give rise to "measured" 
upper limits that can be substantially lower than the "sensitiv- 
ity" or "expected upper limit" of the corresponding search. 
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Figure 7: Scatter plot of the reconstructed variables AE 
and Mlfvd for the final event selection obtained in the 
most recent Belle search for r — > /ry. Dots are data events 
and yellow squares represent MC signal events. The inner 
green ellipse is the 2a signal region (where a total of 10 
data events are observed, with an estimated efficiency e = 
5.1%), while the outer red ellipse is the 3a blinded region. 



3.1. Search for r — ► Z7 (/ = e, /x) 

The most recent r — > /17 and r — > e7 results were 
submitted for publication by Belle [18] using a data 
sample corresponding to a luminosity C = 535 fb -1 . 
The r — > /J7 (r — > e7) analyses have a 5.1% (3.0%) 
signal efficiency within a 2D elliptical signal region in 
the M LFV d and AE plane (see Fig. [7] and Fig.(8|. 

After performing a 2D unbinned extended maxi- 
mum likelihood fit for the number of signal (s) and 
background (b) events, Belle finds s = — 3.9^3;2 
(-0.14t^) and b = 13.9±|;jj (5.14±|'J?) in the r 

( r ~~ > e 7) search. Toy MC simulations are then 
used to evaluate the probability of obtaining such re- 
sults 7 and to estimate the 90% CL upper limits on 
the number of signal events: Ng L (r — > /X7) = 2.0 and 
Ng L (r — > ej) = 3.3. The branching ratio upper lim- 
its following from these results amount respectively to 
B(t /i7) < 4.5 10~ 8 and B(r ej) < 1.2 10~ 7 at 
90% CL. 

The corresponding limits previously obtained by 
BABAR on C = 211 fb -1 amount to B(r -> ^7) < 
6.8 10- 8 [15 and S(r e 7 ) < 1.1 10" 7 [19] at 90% 



7 Belle evaluates that P{s < -3.9) = 25% (P(s < -0.14) = 
48%) for null true signal in the r — > 7x7 (r — > e7) search. 




Mey (GeV/c 2 ) 

Figure 8: Scatter plot of the reconstructed variables Ai? 
and Mlfvd for the final event selection obtained in the 
most recent Belle search for r — > e7- Dots are data events 
and yellow squares represent MC signal events. The in- 
ner red ellipse is the 2a signal region (where a total of 
5 data events are observed, with an estimated efficiency 
e = 3.0%), while the outer green ellipse is the 3a blinded 
region. 

CL. 

Fig. [9] summarizes and compares the BABAR and 
Belle results discussed above, and the corresponding 
analyzed luminosities. 

3.2. Search for r /^° 5 Irj, lrf 9 IK° S 

Both BABAR [20] and Belle [21] have recently pub- 
lished new results on LFV r decays involving the 7r°, 
rj and 7/ pseudoscalars: r — > Z7r°, Z77, Z^ 7 , where Z is 
separately identified as either an electron or a muon. 
The luminosity analyzed by BABAR is C = 339 fb -1 , 
while Belle used C = 401 fb -1 . In the modes with an 
n meson in the final state both the 77 — > 77 and the 
r] — > 7T + 7r - 7r° decays are used. In the r — > Z?] 7 anal- 
yses the 7/ — > 7r + 7r - ^ (77 77) and T7 7 — > p°7 decay 
modes were included. Fig. [To] shows the distribution 
of the finally selected data events in the Mlfvd and 
AE plane for each of these channels, together with 
the 2a signal box used in the BABAR search. 

The BABAR expected background per channel is be- 
tween 0.1 and 0.3 events. Summing over all ten modes, 
the total expected background within the signal re- 
gions amounts to 3.1 events, whereas 2 events in total 
were observed. 

Belle had previously published a search for an addi- 
tional mode with a neutral pseudoscalar meson in the 
final state (r -> IK% [22 , C = 281 fb" 1 ). 

Fig. [TT] summarizes and compares the BABAR and 
Belle results on all searches for r LFV decays to a 
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Figure 9: Summary and comparison of the current upper limits (dots: left vertical axis) set by BABAR (red) and Belle 
(blue) on B(r — > /jj) and B(r — > ej). The corresponding analyzed luminosities are shown as horizontal lines on a scale 
displayed on the right vertical axis. 
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Figure 10: Selected data (dots) and 68% of signal MC events (shaded region) in the Mlfvd and AE plane for the 10 
LFV decay channels r — > Itt°, In, In' searched for by BABAR; for each channel, the corresponing 2a signal box is displayed. 



lepton and a neutral pseudoscalar meson. 



of decays r -> 3/ [25l [26] 8 , while Fig. [14] shows the 
results published by Belle on several r decays to a 
lepton and a neutral vector meson [24 . 



3.3. Other searches for r LFV decays 



Both BaBar and Belle have searched for other r 
LFV decays. The corresponding results are compared 
in Fig. 12 for the class of decays r ^_lhh! (where h 
and h! = tt z 



Fig. 13 compares the 



results obtained by the two experiments for the class 



8 BABAR has very recently submitted for publication |27| an 
updated search for r — > 31, using a much larger luminosity 
(C = 376 fb -1 ) than previously published by each of the two 
experiments (see Fig. |13| . Upper limits on the branching frac- 
tions are set in the range (4 — 8) 10 -8 at 90% CL. These results 
were not public when FPCP2007 took place, and were not pre- 
sented at the Conference. 
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Figure 11: Summary and comparison of the current upper limits (dots: left vertical axis) set by BABAR (red) and Belle 
(blue) on the branching ratios for r decays to a lepton and a neutral pseudoscalar meson. The corresponding analyzed 
luminosities are shown as horizontal lines on a scale displayed on the right vertical axis. 
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Figure 12: Summary and comparison of the current upper limits (dots: left vertical axis) set by BABAR [23] (red) and 
Belle 24 (blue) on the branching ratios for r decays to a lepton and a pair of charged pseudoscalar mesons. The 
corresponding analyzed luminosities are shown as horizontal lines on a scale displayed on the right vertical axis. 



4. Future Prospects for LFV searches in 
r decays 

The estimated physics reach of the full data sam- 
ple that will be collected by BaBar and Belle till the 



end of their foreseen running, based on projections 
from existing analyses, depends on the residual back- 
ground level. We express the experimental reach in 
terms of the "expected 90% CL upper limit" and, for 
brevity's sake, refer to this as the "sensitivity" . In the 
absence of signal, for large A/"bkd? ^9o L ~ l-647Vbkd> 
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Figure 13: Summary and comparison of the current published upper limits (dots: left vertical axis) set by BABAR [25] 
(red) and Belle [26] (blue) on the branching ratios for r decays to three charged leptons. The corresponding analyzed 
luminosities are shown as horizontal lines on a scale displayed on the right vertical axis. 
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Figure 14: Summary of the current upper limits (dots: left vertical axis) set by Belle [24] on the branching ratios for r 
decays to a lepton and a neutral vector meson. The corresponding analyzed luminosity is shown as a horizontal line on 
the scale displayed on the right vertical axis. 



whereas for small A^bkd a value for is obtained 

using the method described in [28J. So, for TVbkd ~ 0? 
^90 L ~ 2.4. Reducing the background below a hand- 
ful of events doesn't greatly improve the expected 
limit if significant efficiency is lost in the process, 
which is why it is common to see experiments report- 
ing the expected backgrounds to be small (i.e. a few 
events), but rarely below 0.1 of an event. 



A worst-case scenario is obtained if identical analy- 
ses to those already published by BaBar and Belle are 
repeated, as is, on the increased data sample: in that 
case the expectations then simply scale as y^bkd/A 
which, for large iVbkd? scales as 1/y/C. A best case 
scenario would take the current expected limit and 
scale it linearly with the luminosity. This is equiv- 
alent to a statement that analyses can be developed 
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maintaining the same efficiency and backgrounds as 
the current analyses. 

For r — > Z7, there is an irreducible background from 
r — > Zz/i/ for e + e~ — > r + r~ events accompanied by 
initial state radiation (ISR). The ISR photon can be 
combined with the lepton I to form a r — > Z7 candi- 
date that accidentally falls in the signal region in the 
^lfvd — plane. Such a background is irreducible 
in the sense that it arises from an e + e~ — » r + r~ pro- 
cess with a well measured lepton and 7 in one of the 
r hemispheres. In the existing BABAR analyses, these 
events account for approximately one fifth of the to- 
tal residual background. Scaling with this irreducible 
background only, one would expect an upper limit for 
B(t — > Z7) which ranges between 1 and 2 10 -8 from a 
complete combined BABAR and Belle data set. 

The situation for the other LFV decays, r — > 31 and 
r — > is even more promising, since these modes 
do not suffer from the aforementioned backgrounds 
from ISR. In this case, one can project sensitivities 
assuming iVbkd comparable to backgrounds in exist- 
ing analyses for approximately the same efficiencies. 
These yield expected limits at the 10 -8 level with the 
complete BABAR and Belle data set. 




Figure 15: BABAR and Belle current experimental limits 
on B(t — > /jlj) compared to a range of theoretical predic- 
tions for this branching ratio in a SUSY $0(10) scenario 
of BSM physics [9 . The best experimental limit available 
before the two asymmetric B-factories started to play a 
role is shown by the uppermost dark green horizontal line 
("PDG2000 limit"). 



5. Conclusions 

Searches for LFV in r decays are an optimal hunt- 
ing ground for BSM physics, which is complementary 
to possible LHC discoveries: observation or non ob- 
servation of LFV processes in the charged sector can 
significantly constrain theory parameter space. 

BABAR and Belle have looked for signals of LVF in 
many exclusive r decay modes and, although no evi- 
dence of LFV has yet been found, limits have pushed 
into the 10 -8 region where parameter space in some 
BSM scenarios is already being constrained. An ex- 
ample of the constraining power of the current exper- 
imental limits is shown in Fig. [l5j where the BABAR 
and Belle limits on B(r — > /ij) are superimposed to 
a scatter plot of theoretical predictions obtained in 
a SUSY 50(10) scenario of BSM physics [9 . Also 
shown for comparison is the best experimental limit 
available in the pre-asymmetric-5-factory era. 

It should finally be mentioned that there exist pro- 
posals [29 for Super ^-factories which would generate 
up to a 100 fold increase in the size of the r sam- 
ple compared to those expected from the existing B- 
factories. If such a facility is built, one will be probing 
LFV r decays at the O(10- Q lO" 10 ) level. 
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